Geometric Phase Based Quantum Computation Applied to an NP-Complete Problem 
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We present a new approach to quantum computation involving the geometric phase. In this 
approach, an entire computation is performed by adiabatically evolving a suitably chosen quantum 
system in a closed circuit in parameter space. The problem solved is the determination of the 
solubility of a 3-SAT Boolean Satisfiability problem. The problem of non-adiabatic transitions to 
higher levels is addressed in several ways. The avoided level crossings are well defined and the 
interpolation can be slowed in this region, the Hamiltonian can be scaled with problem dimension 
resulting in a constant gap size and location, and the prescription here is sufficiently general as 
to allow for other suitably chosen Hamiltonians. Finally, we show that with n applications of this 
approach, the geometric phase based quantum computation method may be used to find the solution 
to the 3-SAT problem in n variables, a member of the NP-complete complexity class. 

PACS numbers: 03.67.Lx, 03.65. Vf 



I. INTRODUCTION 

Geometric phases have been proposed in the context 
of quantum computation, specifically in the proposal of 
fault-tolerant gatesQQQ- Here we discuss the possibil- 
ity of performing an entire computation using a single 
application of the geometric phase. This approach uses 
the geometric phase phenomenon in conjunction with the 
adiabatic approach of Farhi, et. al. j4|- With a single 
application, the problem solved is determining the sol- 
ubility of a Boolean satisfiability 3-SAT problem in n 
variables. With n applications of this approach, the ge- 
ometric phase based quantum computation method may 
be used to find the solution to the 3-SAT problem it- 
self. This last problem is of the NP-complete complexity 
class 0. 



II. ADIABATIC QUANTUM COMPUTATION 

Farhi and coworkers presented an approach to quan- 
tum computation, Adiabatic Quantum Computation 
(AQC). The AQC approach involves an adiabatic in- 
terpolation from the ground state of an initial Hamil- 
tonian (Hi) to the ground state of a final Hamiltonian 
(Hp) that encodes the solution to the problem of interest. 
The method involves adiabatically evolving the Hamilto- 
nian H{z) = zHi -|- (1 — z)Hp, while remaining in the 
ground state. In this approach, the problem solved is an 
NP-Complete problem, the 3-SAT Boolean Satisfiability 
problem containing m clauses in n variables. In our con- 
struction, each eigenstate of the N = 2"" dimensional Hp 
encodes a possible solution to the problem. The energy 
of each level of Hp corresponds to the number of clauses 
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that have been violated by that possible solution, with 
the ground state E=0 level corresponding to a solution 
of the problem (0 violated clauses). 

In the context of AQC, it was shown that in the case of 
a computationally hard instance of the 3-SAT problem, 
there exists a parameter region in which the spectrum 
is irregular, the spectral fluctuation distribution taking 
the form of the Wigner distribution!^. In these regions, 
it was shown the interpolation was susceptible to non- 
adiabatic transitions from the ground state leading to an 
interpolation time that scales exponentially with prob- 
lem size. The question arises whether similar AQC ideas 
could be implemented without interpolating through the 
irregular spectral region. Namely, by adding a second 
parameter dependence and interpolating around an ap- 
propriately modified final Hamiltonian, the geometric 
phase may be used and a computation performed without 
traversing through a problematic region. 



III. GEOMETRIC PHASE 

BerryQ noted that an eigenstate having a parameter 
dependence acquires, in addition to the familiar dynam- 
ical phase, a geometric phase (Berry phase) when the 
state is adiabatically evolved in a cyclic manner near a 
level degeneracy, termed a diabolical point. In the spe- 
cial case of a real Hamiltonian, cyclic evolution of the 
parameters produces a phase of tt, an overall sign change 
in the eigenstate being transported. 

In the current approach, we use a construction similar 
to that of AQC - although instead of interpolating in one 
parameter (z) from Hi (z = 1) to Hp [z = 0), we use 
two parameters (x, z) and evolve the system in a closed 
circuit in the two dimensional parameter space about the 
point (0,0). To illustrate the method, we show how the 
Berry phase phenomena yields a test for solubility of the 
3-SAT using Hp itself. Hp encodes a particular 3-SAT 
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problem instance in n variables and is constructed easily 
without a priori knowledge of the solution^. 

The Haniiltonian for this problem is constructed in 
such a way that a level degeneracy only exists at the 
point (0,0) if a solution exists, and no level degeneracy 
exists if a solution does not exist. When the system is 
transported along the closed path encircling the point 
(0,0), the state accrues a Berry phase if a solution exists 
and accrues no phase if a solution does not exist. Since 
the Hamiltonian is real, the accrued phase becomes tt in 
the event of the existence of a solution^. 



IV. GEOMETRIC PHASE BASED QUANTUM 
COMPUTATION 

The Hamiltonian we use for this problem is the direct 
sum of two Hamiltonians with the addition of interac- 
tion terms. The first Hamiltonian, Ha, is constructed by 
adding |I to the Hp above: Ha ~ Hp + |I, which has 
N = 2" dimensions, where n is the number of variables 
in the 3-SAT problem instance. We construct Ha in such 
a way that relies on no a priori knowledge of which state 
of Hp is the solution state. This definition of Ha is se- 
lected for certain simplifications, although many forms 
for Ha may be chosen. In fact, it is worth noting that 
Ha can be chosen to be Hp identically. In the follow- 
ing, we demonstrate the method for a soluble 3-SAT in- 
stance, with the generalization to non-soluble cases being 
straightforward. We arrange the basis states of Hp such 



that i/p,, 7^ 0, 1 < i < iV-1; and Hp 



(the solution 



state). The second Hamiltonian under consideration, in 
addition to Ha, is the Hamiltonian Hb = — f , having di- 
mension 1. We define the unperturbed Hamiltonian for 
the problem to be If^"^ = Ha^Hi,. The unperturbed 
states lipi) jlipa) , and \ipb) are defined by: 



I i^a) 

h'^"^ I A) 



= H 



= H 



(0) 



\A) = i^ + Hp^,)\^P.) ,it<N) 



I V'a) = J I ^a) 



(1) 



The final form for the Hamiltonian used in the geometric 
phase computation includes coupling between Hj, and all 
states of Ha through the additional parameter x. In the 
unperturbed basis of Eq. |^ , the general matrix represen- 
tation of H that determines the 3-SAT solubility through 
the geometric phase phenomenon is 



Hi, 



Hau 



H„ 



[i < N) 



Hnn — Hgnn — — + H, 



'-PNN 



HiN = Hn, = x {i < N) 
Hn+i,n+i = = , 



(2) 



where no particular problem instance has yet been chosen 
for Hr,. 
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FIG. 1: Spectrum generated by Eq. (1) for varying z: (a),(b); 
varying x: (e),(f); the level crossing: (c),(d); and avoided level 
crossings: (g). Note: scale varies between rows. Note also 
(g),(h) are identical to (e),(f) with a change of scale. The 
path in parameter space traverses frames (f) (r to 1) —* (a) (r 
to 1) ^ (e) (1 to r) ^ (b) (I to r) ^ (f). 



The problem to be solved is determining the existence 
of a solution to the 3-SAT problem, or 3-SAT solubility. 
Evolving the ground state of Eq.© adiabatically in a 
closed circuit in the parameters (x,z) results in one of 
two possible outcomes: no geometric phase in the case of 
no solution, or a geometric phase = tt in the case of the 
existence of a solution to the 3-SAT problem. 

Our numerical analysis uses the specific case of n=7, 
N=128. In addition, we choose a specific problem in- 
stance encoded by Hp. We choose the highly degenerate 
case of Hp having a solution state (E=0) and N-1 degen- 
erate states (E=l). Specifically, Hp^. = 1, i < N and 
HpNN — 0- Although somewhat unrealistic from a prob- 
lem standpoint, this case produces the smallest minimum 
gap size at the avoided level crossings (see Fig. Ig) where 
a non-adiabatic transition may occur, and so represents 
the worst case situation. This is discussed further below. 

We start the system in the ground state of the sys- 
tem at (x,z)=(0,l) (see Fig. Ih); this corresponds to 
the unperturbed state \ipb), with the eigenvalue equation 
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H{0, 1) li^biO, 1)) = - J \M0, !))• This state is then adi- 
abatically evolved around the following parametric cir- 
cuit: {x,z) = (0,1) ^ (-1,1) -> (-1,-1) ^ (1,-1) ^ 
(1,1)^(0,1). 

The approach may be viewed graphically in the follow- 
ing manner. The spectrum generated by Eq.jSJ with the 
particular instance of Hp so defined is shown in Fig. 1. 
This choice of Hp produces a spectrum of three energies 
levels. The highest energy sheet has a degeneracy (N-1), 
and the two lowest are nondegenerate. These lowest two 
energy sheets are responsible for the generation of the 
geometric phase. 

If the 3-SAT problem encoded in Hp has a solution, 
then the ground state energy of Hp is and there will 
be a degeneracy in the two lowest energy sheets at the 
origin in parameter space, {x = 0, z = 0)(see Figs. Ic, 
Id). When the ground state is evolved adiabatically in 
a closed circuit in parameter space about this point, this 
will give a Berry phase = tt (see Figs, la, lb, le, If). 
If Hp is chosen such that no solution exists, there is no 
degeneracy, and the closed circuit does not encircle a di- 
abolical point resulting in no Berry phase. 



(2) (2) 

by Aa and AJ^ ' . The second order correction is given by 
the famihar formula [10| A^^ = J2n^rn iWTiW' where 
H' is the perturbation matrix containing only off diago- 
nal terms equal to x {H- ~ ^tv-i-i i = * ^ We 

(2) 

begin with A^ . The sum in this ease can be split into 
two terms, one involving the (N-1) states with energy 

^ and one involving the state with energy Eji^\ In 
both cases the matrix element yields | i/^ p= x'^{n < 

N). Ai^-* simplifies directly since H'^ ^ = for z 7^ iV + 1, 
and the sum reduces to one term. Evalution of these 
sums, using Eq.lO, yields 



A(2) = -2x^ 
(2) _ -2(jV-l)x^ 2x^ 
" ^ 2 + z 



(4) 



Inserting N=128, z=-l into Eqs.®,® the energies to 
2nd order are: 

E„ = ---2x^ 



AVOIDED LEVEL CROSSINGS 



Eh = - - 252x^ 



(5) 



We now consider the problem of avoided level cross- 
ings: gap size and location. 

In order to encircle the origin and accrue the Berry 
phase, the state must remain on the ground state en- 
ergy sheet and not make a transition at the avoided level 
crossing (eg. by way of Landau-Zener transition [3)(see 
Fig. Ig). Simulations have shown that the avoided level 
crossing gap size produced by Eq.Q with Hp defined 
above decreases exponentially with increasing problem 
size. This is similar to difficulties in the AQC approach 
1^. Even so, unlike other approaches, the location of the 
avoided level crossings is well defined, can be predicted to 
some extent, and the interpolation slowed in that region 
to prevent a level transition. 

We gain insight into the Hamiltonian in Eq. jSJ by cal- 
culating the approximate location of the avoided level 
crossings using perturbation theory. We consider the 
Hamiltonian in Eq.Q with Hp defined above evaluated 
at z=-l and x small (corresponding to Fig. Ig). In this 
case, the unperturbed states are given by Eq.QJ, taking 
X as the perturbation parameter. As x does not appear as 
a diagonal entry in the perturbation matrix, there is no 
first order correction. The second order correction is eas- 
ily evaluated using the following unperturbed energies: 



E. 



(0) 



H 



(0) 



Hh = 1 



{i < N) 



E, 



(0) 



H 



N+1,N + 1 



We are interested in the second order perturbative cor- 
rections to E^^ and Ej^^\ and we denote these energies 



This is consistent with Fig. le and Fig. Ig. The differ- 
ences are attributed to higher order perturbations,which 
become more pronounced with increasing x. The in- 
tersection of these two parabolas gives the approximate 
location of the avoided level crossings to second order, 
Xgap « ±0.045 (see Fig. Ig). 

Earlier we chose the Hamiltonian Hp such that Hp.. — 
1, i < N and Hp^,^ = 0. The reason results from the 
perturbation theory calculation above: each level above 
contributes to the level repulsion on state \il)b) (see Fig. 
Ig). In this case for | x |< Xgap, (N-1) levels contribute 
to the level repulsion downward, and one state (iV'a)) 
repels the state upward. In the general case, the sum is 
dominated by the closest levels. The Hp chosen has N-1 
levels located at the closest energy corresponding to a 
non-solution (E=l). This results in the situation having 
the strongest level repulsion, and in this sense represents 
the worst case situation in terms of the avoided level 
crossings. Other choices of 3-SAT problem instances for 
Hp result in weaker level repulsion. 



VI. SCALING 

In addition to varying the interpolation speed near the 
avoided level crossings, the problem can be scaled with N 
such that the minimum gap size remains N independent 
for all problem sizes. We choose the scaling such that 

('2) 

A}j is N independent for large N in Eq. I^J . This may be 



accomplished by scaling z or by scaling x appropriately. 
If z is chosen to scale with N such that z Nz, the per- 

approach a constant 



f2l 12 

tubative corrections Aa and A^ 
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FIG. 2: Energy level plots showing the crossing and avoided 
level crossing along the path in parameter space including the 
scaling in Eq. Not shown are the (N-1) degenerate levels 
at Energy=z/4 + A^Hp,, « 128. 



for large N. Explicitly, the Hamiltonian containing this 
scaling is the following (cf. Eq.©): 

Ha = Ha,, = z/4 + NHp^^ 1 <i< N 
Hnn = Ha^N = 2/4 + NHpj^j^ 

HiN = Hn, ^ X 1 <i < N 

Hkk = Hb = ~z/4: k = N + l. (6) 

Wc have used n=7, N=128 in the numerical simulation 
of Eq. ©, with the results shown in Fig. 2 (cf. Fig. 1). 
The scaling was checked numerically for 3 < n < 9. For 
these values, a constant minimum gap size (Ag^p ~ 0.5) 
and constant gap location {xgap ~ 0.0) was found. 

Analogously, we may choose N scaling for x, such that 



(2) i 

Aa and AJ,^'' is approximately constant for large N. This 
may be accomplished by choosing to scale with -ir , or 



.(2) 



VII. SOLVING THE 3-SAT PROBLEM 

The geometric phase based approach above determines 
the existence of a solution of the encoded 3-SAT prob- 
lem in n variables. We now show that upon repeated 
iterations of this approach, the solution to the 3-SAT 
problem may be searched for and the solution found in 
n iterations. 

As described above, if a solution state exists a geomet- 
ric phase will be generated. Initially, each eigenvector 



{(pi) of Ha encodes a possible solution to the 3-SAT prob- 
lem. The collection of these eigenvectors forms the basis 
of Ha and encodes all possible solutions. This represents 
the search space under consideration. From among this 
set we denote the state \(j)a) as the solution state, if it 
exists. 

Using this Ha, wc construct H as in Eq.|(21). If the 
Berry phase is generated, then the solution state exists 
in the search space. If no phase is generated, then no 
solution exists and the search is complete. 

In the case of the existence of a solution, we then divide 
the search space in half, such that Ha = h'^^ ®Hf\ 
with the basis of each subspace forming a solution space, 
each containing ^ solutions. 

Substituting h'^'' for Hp in Eq.Q, we construct a new 
Hamiltonian, H^^\ having y + 1 dimensions. Performing 
the geometric phase approach with this H^^\ the system 
either accrues a phase, in which case the solution state 

exists in the basis of Ha"^ ; or the system accrues no phase, 
implying the solution lies in the basis of Hj^\ 

The Hamiltonian containing the solution state in its 
basis is then subdivided and used in the next iteration in 
a similar manner. For example, if Ha^ generates a phase, 
then the solution exists in the basis of Ha^ which is then 



subdivided into H^^^ = H^^'^H'^'. Substituting H, 
for Hp in Eq.©, wc construct a new Hamiltonian, H'^'^\ 
having ^ + 1 dimensions. A single application of the 
geometric phase approach determines the solution to be 
in the basis of iji^^ or H^p\ 

This iterative procedure is a sorting algorithm which 
halves the search space at each iteration, and locates 
the half containing the solution in a single evaluation 
(the geometric phase based procedure). Starting with N 
states, the search is completed in log2(iV) evaluations. 
In the case of the 3-SAT problem in n variables, the ini- 
tial search space is = 2" yielding the solution in n 
evaluations. 



r(2) 



(2) 



VIII. CONCLUSION 

We have illustrated a general approach to solving 
NP-completc problems using the geometric phase phe- 
nomenon in conjunction with the Adiabatic Quantum 
Computation method. 

We first address the problem of whether a Boolean 
Satisfiability problem (3-SAT) is soluble - if so, the adi- 
abatically evolved eigenstate acquires a geometric phase; 
if the problem is not soluble, the state does not acquire 
a phase. 

The problem of avoided level crossings having small 
gap sizes is addressed in several ways. First, the crossings 
are isolated and well defined, and so the interpolation 
can be slowed in this vicinity. Second, the Hamiltonian 
can be scaled with N (see Eq.®) to make the gap sizes 
independent of problem size. Third, as mentioned earlier. 
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the Hamiltonian (Eq.Q) is not the only choice which 
leads to a successful approach. This is a particularly 
simple construction, but there is a large degree of freedom 
in choosing the Hamiltonian that governs the approach. 
Other Hamiltonians can be found that also adhere to 
additional experimental requirements or parameters of a 
different problem. 

Repeated iteration of this procedure results in a search 
of the set of N basis states encoding the possible solutions 
that completes in log2(A^) iterations. For a problem in- 
stance having n variables, the solution of the 3-SAT can 
be found in n iterations of the geometric phase based 
approach. In this way, we show the geometric phased 



based approach yields the solution to the 3-SAT, an NP- 
complete problem. 
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